When translocation dynamics becomes anomalous 
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Recent single molecule experiments probing the passage process of a short single-stranded DNA 
(ssDNA) through a membrane channel (translocation) allow to measure the passage time distribu- 
tion. Building on a recent modelling approach (D. K. Lubensky and D. R. Nelson, Biophys. J. 77, 
1824 (1999)), which has been demonstrated to be valid for chains of up to ~ 300 nucleotides and 
therefore well applies to the system we have in mind, we discuss the consequences if the associated 
dynamics is not of Markov origin, but if strong memory effects prevail during the translocation. 
Motivation is drawn from recent results indicating that the distribution of translocation times is 
broader than predicted by simple Markovian models based on Brownian motion. 
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The translocation of biomolecules through membrane 
pores (channels) is one of the most vital processes 
within or across biological cells, serving both d eliv- 
ery and signalling purposes j Alberts et ah L 17994). In 
(bio)chemistry, forced translocation is used in selec- 
tion/purification of larger molecules, and in medicine, 
it plays an important role in drug delivery. Whereas 
the translocation of short, inflexible molecules is primar- 
ily determined by the properties of the pore (energy- 
driven transport, sticking events within the pore, etc.) 
and the difference of the chemical potential between the 
cis and trans sides of the pore, semiflexible and flexi- 
ble molecules, in addition, have to cross an entropy bar- 
rier while being (partially) confined within the channel 
llLubenskv and Nelson lll999UMuthukumar L l200lj.ll999t 
ISlonkina and Kolomeiski I . l2003t ISung and Park I . Il996). 
In the presence of a high external bias and for the rather 
short chains used in typical experiments, the entropic 
slowdown as well as the other interactions between chain 
and channel wall become neglig ible, the passage being 
domi nated by the applied drift l)Lubenskv and Nelson I 
1999). In what follows, we develop a scenario accord- 
ing to which the translocation dynamics is governed by 
slowly decaying memory effects, leading to a different 
behaviour in the distribution of passage times which we 
believe can be measured experimentally. 

Experimentally, the translocation of ssDNA can be 
observe d on a single mole c ular level, both v o ltage- 
driven llAkeson et al. L Il999t iKasianowicz et al. L Il996t 
iMeller et al. L l200l[) and in the absence of an exter- 
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nal electric field l)Bates et al. L l20f)3j) . In such single- 
molecule translocation assays, fairly short chains are 
used, with some 60 bases corresponding to approx- 
imately 12 persistence len gths, or 6 Kuhn lengths 
l|Frank-Kamenetskii L Il997t) . The width (~ 50 A) of 
the membrane amounts to about one persistence length 
(~ 40l) of the ssDNA. A good measure for the transloca- 
tion process is the distribution of passage times, i.e., the 
statistics of time spans the chain needs to cross from the 
entry (cis) side to the exit (trans) side of the pore. In the 
results, one observes two (or three) different time scales: 
the shortest corresponds to chains which retract from the 
pore back to the cis side, before completing the passage 
through the pore; the other (one or two) correspond(s) 
to real passage times (if there are two peaks, this can be 
explained by different ori entations of the chain in respe ct 
to the passage direction l)Lubenskv and Nelson I fl999)). 

In a recent experiment l)Bates et al. I l2003j) . it was 
noted that the first passage time distribution contains 
non-negligible contributions over a large time-range even 
in the presence of a low driving voltage, a case in which 
a Markovian model would predict exponentially fast de- 
cay. This may well indicate that additional mechanisms 
so far neglected play a role in the translocation dynam- 
ics, which might effect long-tailed first passage time dis- 
tributions, and therefore imply a possible modelling by 
assuming a non-Markovian behaviour of the system. In 
this note, we construct a framework in the limit of strong 
non-Markovian effects, taking into account anomalous 
translocation dynamics through long-tailed memory ef- 
fects. Giv en the accuracy of the newly reported experi- 
ments in l)Bates et al. 1 12003^ . it might well be possible 
to resolve such effects in log-log analyses of the presently 
available, or future data. We collect a number of possible 
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sources for such anomalous dynamics. 

In the presence of a bias field and for chains with < 300 
nucleotides, the translocation dynamics in the Markov 
limit has been shown to follow the Smoluchowski-type 
equation ijLubenskv and Nelson I [199^) 

dP ( d d 2 \ 

where P{x, t) is the probability density function (pdf ) 
to find the chain at position x at time t, and v 
and K are the associated drift and diffusion constants, 
which may be determined from more microscopic models 
l|Lubenskv and Nelson L Il999|h The translocation pro- 
cess without retraction from the pore can thus be de- 
scribed by the first passage time distribution F(t) from 
the point x = L to x = 0. In the presence of th e external 
drift, this leads to the result [see, for instance, ijRedner I 
2001)] 



Fi(t) 
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Note the sharp exponential decay for longer passage 
times. In this picture of Markov advection-diffusion, 
the mean first passage time from the origin to an ab- 
sorbing boundary a distance L away is given by T\ = 
J™tFx(t)dt = L/v, i.e., the statistical mean corresponds 
exactly to a classical linear motion with the drift velocity 
v. Thus, the influence of the diffusivity in this average 
becomes negligible and the motion can be characterised 
by the mean (x) = vt. 

This picture dramatically changes in the presence of 
long-tailed memory, effected by a waiting time distribu- 
tion 



t l+a> 



(0<a< 1), 
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according to which interruption times of the transport 
process in a multiple trapp ing model are distributed 
l)Metzler and Klafter Ll2000lafl . i.e., the transport process 
is being stalled successively, and the sticking intervals 
follow (J3J). In this situation, the relation for the mean 
gets modified to (x) — v a t a , where v a is an appropri- 
ately generalised velocity, and the corresponding gener- 
alisation of equation (JIJ is the fractional Fokker-Plan ck- 
Smoluchowski equation ijMetzler and Klafter I . l2000lcJ) 
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with K a = Kv a /v, and the fractional Riemann-Liouville 
operator D\~ a P{x,t) = w% Jo 0J^ dt The 
waiting tim e in expression is scaled to the internal 
time scale r ijMetzler and Klafter ll2000|) . A typical feau- 
ture in this anomalous case is that the maximum of the 
pdf due to the strong persistence of the initial condition 
(i.e., the probability J Q ip(t)dt of not moving), remains 




Figure 1 Translocation probability 1 — S(i) for long-tailed 
waiting time wi th index a = 1/2 (full line), in comparison 
to the Markov counterpart (dashed) on a log-lin plot. In the 
inset (lin-lin plot), the much slower increase of 1 — S(t) in the 
non-Markovian result is even more distinct. In the plot, the 
dimensionless quantities, L — 5, and v = 1 were chosen. 



at the initial location x = 0, and the mean first passage 
time diverges: T a = f Q F a (t)tdt — > oo. The associated 
first passage time distribution in the presence of can 
be determined from the classical result Recalling 
the scaling relation for the pdf P(x, u) in Laplace space 
(P(x, u) = L°° P(x,t) exp(-ut)dt) betwee n the solutions 
of and (QJ (|Metzler and Klafter 1 12000|) . the following 
scaling holds between the Markov survival probability 
Si(t) = J Q F(t)dt and its non-Markov analogue, 



S a (u) =u a - 1 S 1 (u a ) 



(5) 



in rescaled variables. This scaling relation can be rewrit- 
ten in terms of the generalised Laplace transformation 



/•OO 

S a (t)= / E s (s,t)S 1 {s)ds, 
Jo 



(6) 



where the Laplace transform of the kernel E s (s,u) is 
given by the one-sided Levy dis tribution E s (s,u) = 
u"" 1 exp (— su a ), and therefore ijMetzler and Klafter I 
2000) 



E s (s,t) = - 



(-iy 
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In particular, for a = 1/2, one obtains E s (s,t) — 
(jtt)~~ x l 2 exp (— s 2 /(4i)). The pronounced difference be- 
tween the Markov result and its non-Markov analogue in 
the presence of the long-tailed distribution for the ex- 
perimentally measured quantity 1 — S(t), i.e., the translo- 
cation probability, is displayed in figure ^ 

It is straightforward to show that there exists a similar 
scaling relation for the first passage time density: 



F a (u) =F!(u a ). 



(8) 



From this relation, and the small-it expansion F-\ (u) ~ 
exp (— Lu/[2v]), one can by Tauberian theorems ijFeller I 
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Figure 2 First passage time density F a (t) for a = 1/2 (full 
line; v = 10, upper curve, and v = 5) and Markov case 
(dashed, « = 10). Not e the higher peak in the Markov 
case, indicating that the tail falls off faster than in the case 
with long-tailed waiting times. The inset shows the double- 
logarithmic plot with the —3/2 power-law asymptotics, for 
v = 10, 5, and 2.5. 



1968) deduce the long-time behaviour of the first passage 
time density, 

FoW-r 1 -". (9) 

Thus, in contrast to the Markov case, in the presence of 
long-tailed waiting times, the density F a itself exhibits 
a long-tail, and the mean first passage time diverges. In 
this case, in analogy to equations and ©, the trans- 
formation from Fi to F a , F a (u) = / °° exp(— u a t)Fi(t)dt, 
can be rewritten as F a (t) = E F (s,t)Fi(s), with 

VM-it-g^-d*)"- < 10 » 

t Till (—cm) \t a J 

For a = 1/2, one infers E F {s,t) = s (47rf 3 )~ 1/2 
exp (— s 2 /(4i)) } In figure [3 we depict the functional 
behaviour of the first passage time distribution for an 
external bias, in comparison to the Markov case, reveal- 
ing the distinct inverse power-law tails in the presence of 
©■ 

Conversely, in the absence of a drift, the reflect- 
ing boundary at x = cannot be a priori ne- 
glected (the retraction of the chain towards the cis 
side is no longer suppresse d by the bias, compare 
l|Lubenskv and Nelson L Il999|) ). and the corresponding 
first passage time distribution follows from the eigen- 
value problem, such that the survival probability in the 



anomalous case is given by a sum over Mittag-Leffler 
decays- S(t) = § E~=-oo ("^j^") , 

which in the Markov limit a = 1 reduce to the exponen- 
tial function, and therefore a considerably faster decay 
l|Metzler and Klafter I . l2000b(L In particular, the long- 
time behaviour of the first passage distribution F(t) ~ 
ensues, i.e., the asymptotic behaviour is analogous 
to the one derived from the case with drift (in contrast to 
the drift-free case without the reflecting barrier, in which 
F(t) ~ t- a ' 2 - r ). 

In both cases with and without drift, the first passage 
time distributions in the generalised case with long-tailed 
memory based on the waiting time form (0, the cor- 
responding F(t) exhibits a distinct maximum, its func- 
tional form being rather similar to the corresponding 
Markov limit, when viewed on a linear scale. On a 
double-logarithmic scale, the power-law F(t) ~ t _1_Q 
differs from the much faster, exponential decay in the 
Markov limit. In particular, for both cases drift-free and 
with drift, the same power-law behaviour is expected, an 
additional consistency check in experiments. The latter 
statement, however, only holds if the mechanism effecting 
the broad form J3J), and the value of a in particular, is not 
affected by the magnitude of the exte rnal field. There are 
indications from the recent studies l|Bates et al. 1 l2003|) 
that this may actually be the case: if the external drift 
is increased, the stalling events become less pronounced. 
This can influence a, but it could also introduce a cut- 
off at some time t c in the waiting time distribution tp(t), 
and therefore cause a system response which is closer to 
the Markov case than in the absence of the drift. We 
stress that in the presence of long-tailed first passage 
time distributions, the most probable passage time cor- 
responding to the value at which F(t) peaks, becomes a 
rather meaningless quantity, as the mean first passage 
time diverges . We also note that the fit in reference 
((Bates et al. I [2003) of the passage time distribution by 
two exponentials with significantly different time scales 
in our approach becomes replaced by a continuum distri- 
bution of relaxation times. 

Let us now list a number of potential sources for the 
waiting time distribution ?/>(£): 

(i) During the passage, sticking events caused by pore- 
chain interactions may occur. As these events would be 
expected to be correlated with the co-operative motion of 
the pore molecules, and possibly the translocating chain 
itself, a good guess would be that this causes power-laws 
of the form lEH ), e.g. as observed in ligand rebinding in 
proteins l)Glockle and Nonnenmacher I Il995l) . 

(ii) It is known from patch clamp measurements on 



1 Note that in the case discussed here no turnover between two 
power-laws, from t a ~ 1 to t~ a_1 results, in contrast to the prop- 
erties of the cases discussed in H. Scher and E. W. Montroll, 
Phys. Rev. B 12, 2455 (1975); H. Scher, G. Margolin, R. Met- 
zler, J. Klafter, and B. Berkowitz, Geophys. Res. Lett. 29, 
10.1029/2001GL014123 (2002). 



2 The Mittag-Leffler function E a (z) = Eo° z n /V(l + an) is the 
'natural generalisation' of the exponential function e z = E±(z). 
E a (—ct a ) interpolates between an initi al stretched exponen- 
tial and a final inverse power-law ~ t~ a iMctzlcr and Klafter , 
2000). 
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single fast chloride channels ijBlatz and Maglebv I Il986() , 
that ion channels open and close repeatedly. For the dis- 
tribution g(t) of duration times of such channel gating 
events, it was typically found that it follows a modu- 
lated power-law trend, g(t) = Aft)^ 13 , j3 > 0, where 
the modulation factor A can follow logarithmic oscilla- 
tions before eventually being cutoff by an exponential 
llNonnen macher and No nnenmacher I Il989j) . or be con- 
stant ijMillhauser et al. Ul988|) . Within a finite time win- 
dow, both are indistinguishable. It is therefore fair to say 
that gating events in a given time window in single ion 
channels follow power-law statistics, and typical values 
for (3 are around 1.6. The distribution g translates into 
our waiting time distribution ip(t) from equation J2J with 
a = /3-1, 

(hi) For longer ch ains, Chuang et al. argued 
l|Chuang et al. I . |200ll) . that the diffusion of the chain 
becomes anomalous. Naively viewing the translocation 
as a waiting time process during which the monomers in 
the pore channel have to wait until they are given way by 
the vicinal monomers, and so on, creating a non-Markov 
process which, on some coarse-grained l evel, may well b e 
described by equation @ , compare also l|Douglas 1 (2000). 

This list of scenarios is not meant to be complete. 
However, one might suspect that the sticking scenario 
(i) is most liable to be affected by the strength of the 
external bias, producing an effect simi l ar to the recent 
experiments reported by ijBates et, al. 1 . 12003). in which 
the dynamics exhibits the abovementioned turnover from 
broad to Brownian motion-type statistics on increase of 
the external bias field. 

In some translocation experiments, apart from the 
sharp initial peak in the first passage time density stem- 
ming from immediately retracting chains back to the cis 
side there occurs another hump similar to the one of the 
translocated chains discussed above. It has been argued 
that this is due to the existence of second characteristic 
passage time, depending on the orientation of the chain 
to the membrane channel i n respect to the cis-trans direc - 
tion ('head or tail first') ijLubenskv and Nelson I Fl999). 
The same effect is expected in the case with long-tailed 
statistics following equation (JjJJ. However, it might well 
be that the associated power-law exponent a is different 
for the two orientations, as the nature of the effective 
interactions giving rise to the long-tailed waiting times 
depend on this head-tail difference. 

One might speculate about the biological relevance of 
anomalous translocation dynamics. On the one hand, 
it might be the outcome of a tradeoff between lack of 
specificity, if the passage is too free and a large variety of 
molecules could pass the membrane, and too high sup- 
pression, which would require active transport through 
the pore, implying a fairly large energy cost for long 
molecules. On the other hand, it might be advantegeous 
to have a large variation in the arrival times of translo- 
cated molecules on the trans side (and thereby very ef- 
ficient retention of un-translocated molecules on the cis 
side). 



We have discussed possible changes arising in the dis- 
tribution of first passage times in biopolymer transloca- 
tion through a membrane channel, and listed a number 
of reasons which might give rise to such anomalous be- 
haviour. It should be possible to determine the quantity 
F(t) from experiments to sufficient accuracy, in order to 
be able to distinguish the normal (Brownian) dynamics 
result from its anomalous counterpart in both the pres- 
ence and absence of an external drift. The large qualita- 
tive difference between exponential and power-law forms 
should be easily discernible on a double-logarithmic scale. 
It should, however, be stressed that the onset of the 
power-law trend depends on the strength of the drift, and 
might occur for fairly large times if the drift is weak. We 
finally mention that the proposed long-tailed effects may 
also pertain in other systems, like during the ejection of 
the DNA of bacteriophages from the c apsid through a 
long pipe-like channel into the host cell ijAlberts et al. I 
11994 iMuthukumar L l200l|) . 
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